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Abstract. Resonant soft x-ray Bragg diffraction at the Dy M4^5 edges has been 
used to study Dy multipoles in the combined magnetic and orbitally ordered 
phase of DyB2C2. The analysis incorporates both the intra-atomic magnetic and 
quadrupolar interactions between the 3d core and 4f valence shells. Additionally, 
we introduce to the formalism the interference of magnetic and nonmagnetic 
oscillators. This allows a determination of the higher order multipole moments of 
rank 1 (dipole) to 6 (hexacontatetrapole). The strength of the Dy 4f multipole 
moments have been estimated at being between 7 and 78% of the quadrupolar 
moment. 



1. Introduction 

Resonant x-ray Bragg diffraction (RXD), whicli is Bragg diffraction wliere tlie energy 
of the incident x-rays are tuned to a relevant absorption edge, fias established itself 
as a powerful tool in modern solid state physics to investigate magnetic, orbital and 
charge ordering phenomena associated with electronic degrees of freedom. The data 
obtained from an RXD experiment can be interpreted using a model based on purely 
atomic quantities, neglecting e.g. band and hybridization effects. Models based on 
atomic quantities have the advantage of being suited to the analysis of both x-ray 
and neutron diffraction experiments, as well as absorption, NMR, EPR, muon and 
Mossbauer spectroscopies [T]. The localized nature of rare earth 4/ electrons makes 
them particularly suited to such analysis which conveniently relies upon quantities 
such as atomic tensors and operator equivalents with the many simplifications that 
these bring. In particular, any atomic observable may be decomposed into its {21 + 1)^ 
multipole operator components (with a; = 0, 21 and q — —x, x, where I is the 
orbital angular momentum quantum number corresponding to the electronic shell 
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in question), and any atomic state may be represented by the set of nonvanishing 
multipole components (T^) 0]. (...) denotes the expectation value of the enclosed 
operator on the ground state of the electrons. Multipoles (T^) are parity even. 
Additionally, they have the definite time signature (—1)^, and we refer to time- 
odd multipoles as magnetic and time-even multipoles as charge. In Rare Earth 
compounds, the 4/ electrons are not strongly coupled to the lattice and they may 
retain their degrees of freedom down to very low temperatures, and although ordered 
states of magnetic dipoles are relatively common, multipoles of rank 2 (quadrupole) 
or higher may order under favorable conditions (31 . Orbital ordering is usually labeled 
as simply antiferroquadrupolar (AFQ) or ferroquadrupolar (FQ) since higher-order 
multipoles are not well established by observations, despite the fact that for 4/ 
electrons (/ = 3) multipoles up to rank 6 can be present. In materials with a 
phase attributed to such "hidden" ordering, exactly which multipoles dominate the 
transition to the ordered phase remains a rather controversial topic, as exemplified 
by the compounds Cex^iSti-x^e U, Np02 and particularly URu2Si2 [S]. RXD is 
particularly useful in characterizing such orders, and applications of this technique 
have produced both direct and indirect evidence of high-order multipole moments in 
DyB2C2 mis], Ce:,Lai-:,B6 [9], and Np02 [ID]. 

DyB2C2 has the highest known AFQ ordering temperature of any rare earth 
systemTg at 24.7 K and as such it has attracted much attention, as many other 
compounds where AFQ or FQ ordering is postulated order at temperatures of 5 K or 
less [2]. At room temperature, DyB2C2 exhibits the tetragonal structure VA/mbm and 
undergoes a transition with small alternating shifts of pairs of B and C atoms along 
c at Tq |T1] which reduces the symmetry to P42 /mnm [12]. Below Tjyj = 15.3 K, 
antiferromagnetic (AFM) order with four distinct magnetic sublattices is observed, and 
the moment orientations are postulated to arise from the underlying orbital interaction 
[13] . Additionally the crystal evolves a monoclinic lattice strain at T]\^ [14]. While 
the properties of combined AFM/ AFQ phase (Phase IV) of the similar compound 
H0B2C2 are postulated to arise from an octupole order parameter [14], the case for 
DyB2C2 is not as clear cut due to the presence of the aforementioned lattice strain, and 
the broader variation in its elastic properties at the AFM and AFM-I-AFQ transitions. 

Due to the less stringent experimental constraints, most RXD studies of rare earth 
materials using x-rays probe the L2,3 edges via an El dipole transition {2p — >■ Sd) or an 
E2 quadrupole transition {2p — ^ 4/). In most cases the information obtained from such 
a measurement is relatively indirect in the case of an El resonance, and with relatively 
low intensity in the case of an E2 resonance. Soft Resonant X-ray Bragg Diffraction 
(SRXD) at the Rare Earth M4^5 edges probes the empty 4/ valence states directly 
by exciting a dipole resonance {Zd — > 4/). SRXD experiments at the Dy M4_5 edges 
have shown that the Coulomb (intra-atomic quadrupole) interaction between the 3d 
and 4f shells is significant [15] . The ordered quadrupole moment of the 4f shell in the 
intermediate state reshapes the observed 3d5/2 and 3d3/2 core hole charge density and 
leads to a splitting of the Sd energy levels. This splitting causes interference between 
different pathways of the scattering amplitude, allowing the observation of high-order 
multipole moments of the 4f shell that have rank 3, 4, 5 and rank 6: respectively the 
octupole, hexadecapole, triakontadipole, and the hexacontatetrapole. Quantifying 
these higher order multipoles allows the ground state of dysprosium to be expressed 
in terms of the associated state multipoles, and accordingly a better understanding of 
its electronic structure can be obtained. 

In the present work, we probe the combined AFM / AFQ phase of DyB2C2 using 
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Figure 1. (a) and (b) Diffracted intensity of the (00^) reflection in DyB2C2 
taken with incident a (vertical) and it (horizontal) polarization at 18K and 12K 
respectively, (c) and (d) diffracted intensity of (a) and (b) corrected for absorption 
effects. Note: All data is normalized to the central peak at the M5 edge in the a 
channel, and the ratio between a and n channels is explicitly maintained. 

the same core hole sphtting formahsm previously applied to the AFQ phase [15] , with 
the structure factor extended to include magnetism. Phenomenological changes are 
introduced to mimic differences between magnetic resonances and nonmagnetic ones. 
The presence of magnetic ordering allows the measurement of odd rank multipoles 
(1,3,5) associated with the magnetic order, and it is observed that their presence 
changes the even rank (T^). 

2. Experimental Results 

A DyB2C2 single crystal was grown by the Czochralski method using an arc- 
furnace with four electrodes and cut with (001) perpendicular to the sample surface. 
Subsequently, it was polished and aligned with 26°, 56° or 86° azimuthal angle. Zero 
degree azimuth corresponds to alignment of the b axis in the scattering plane. The 
orbital ordering (OO^) reflection was recorded (without polarization analysis |16| ) at 
the Dy M4,5 edges of DyB2C2 at the RESOXS end-station of the SIM beam-line 
at the Swiss Light Source. The Dy M4_5 absorption edges were characterized with 
fluorescence yield (FY) and electron yield (EY) at RESOXS and the BL25SU beam 
line at SPring-8, respectively. 

The recorded energy profile of the space group- forbidden (OO^) reflection of 
DyB2C2 in the pure AFQ phase [TS] and the AFM -I- AFQ phases are shown in Figure[T] 
(a) and (b) respectively. The integrated intensity was corrected for absorption, and the 
corrected curves arc shown in Figure[l](c) and (d). We draw particular attention to the 
qualitative difference between the high and low temperature spectra, due to the onset 
of magnetic ordering. In the AFQ phase, the energy profile is largely independent 
of the polarization of the incident x-raysto within a scale factor, which is to be 
expected for a resonance that contains only a single order parameter (in this case, 
the quadrupole). Conversely, in the AFM -I- AFQ phase the structure in the energy 
profile for a and tt polarization of the incident x-rays contains quantitative differences 
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in the ratios of all major features. This is indicative of a combination of magnetic 
and orbital scattering which interfere constructively or destructively as a function of 
energy. In the AFM+AFQ phase the intensity measured in the tt channel is generally 
larger than the intensity measured in the a channel, indicative of a strong magnetic 
contribution to the scattering, which is generally expected to be much larger for the 
M5 edge than for the M4 edge. 

3. Theoretical Framework 

It was demonstrated in the previous study that the unique shape of the energy 
dependence of the (OO^) reflection is well described by the splitting of the 3d core 
states. This splitting results in multiple interfering resonators and consequently adds 
structure to the energy profile in analogy to the unusual energy profile of the resonance 
at the M4 edge in Np02 [IH]- The addition of magnetism to the multipolar coulomb 
interaction in the combined AFM + AFQ phase further lifts the degeneracy of the 
3d shell with respect to the magnetic quantum number M. Here J and M refer 
respectively to the total angular momentum and its projection for the 3d hole state 
and —J<M<J. The intra-atomic quadrupole interaction splits the different 
\M\ levels, and the intra-atomic magnetic interaction lifts the remaining degeneracy 
in ±|M|. The relative oscillator amplitudes are determined by the 4f wavefunction, 
which is characterized in terms of multipoles. 

Following the theory for resonance enhanced scattering of x-rays [H [ITl [181 HSl HQ] 
and its extension to the AFQ phase of DyB2C2 [H] we modify the equations for 
the scattered X-ray amplitude to include contributions from the time-odd multipoles 
allowed in the presence of the magnetic order that develops below Tjv- These 
resonances due to odd-rank multipoles will interfere with the even rank ones either 
constructively or destructively. The X-ray scattering amplitude, the square of which is 
proportional to the X-ray intensity da/dfl, is formed from the product of two tensors 

m 

f = Y^X^F^ = J2{-lfX^QF^ (1) 

K K,Q 

with —K < Q < K. Fq describes the electronic response of the sample and is 
constructed from the polarization vectors of the incident and scattered X-ray beams, 
and is discussed explicitly in appendix A. Fq is expressed as 

F^ = {2K +lf/^Y. ^Q.("' (2) 

where the Wigner D tensor DQ^{a, rotates (which refers to local coordinate 
axes), onto the coordinates of the experimental reference frame used for X^q and 
Fq with Euler angles a,/3, and 7. The Dy site symmetry in the magnetic phase of 
2/m restricts q to q = ±1,±2. K = 2 corresponds to quadrupoles (orbitals), while 
K = \ corresponds to magnetic dipoles, and the K = Q contribution corresponds to 
charge order, which is absent at the given Bragg reflection. describes the atomic 
resonance process which is commonly represented by a harmonic oscillator. In our 
previous work the standard theory has been extended so that is a sum of several 
oscillators created by the splitting of the core state. The amplitude A^{J, M) of each 
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oscillator is labeled by the total angular momentum and magnetic quantum 

number M of the core hole created by the resonance process. 

e«( J, M) = [3M2 - J( J + l)]Qj + gjM ■ Hi„t (4) 

where E is the photon energy, Aj is the difference in energy between the degenerate 
3d J shell and the empty 4f states and S/F j is the lifetime of the intermediate state. 

e'^{J,M) is the energy shift of the core levels due to the intra-atomic quadrupole 
interaction Qj and the intra-atomic magnetic interaction created by the unpaired 4f 
electrons Hi„t. 

There are six resonant oscillators at the M5 edge (M = ±|,±|,±i) and four at 
the M4 edge (M = ±|,±i). The oscillator amplitudes interfere and the branching 
ratio between the two edges is defined as a purely real mixing parameter r j. (J, M) 
is constructed from the structure factor of the chemical unit cell ^E*^: 



r ^ ^ X ^ ^ 

where the symbol in brackets is a 3j symbol [H] and i?^(r, x) are reduced matrix 
elements, with r = 0, 1, . . . , 2 J, x — ji^ — r|, . . . , ji^ -I- r| and g + x and r + x are both 
even integers. The relevant reduced matrix elements R^{r,x) are listed in appendix 
C. 

The Dy site symmetry further dictates that A2 — and A^ — A^^, which 

follows from the properties of 4*^. The quantity ^E*!^ is constructed from the matrix 
elements of the multipole operator combined with spatial phase factors and the 
symmetry components corresponding to the 4(c) site that the Dy occupies, discussed 
in detail in Appendix B and [22] ■ For the (00^) reflection, this is 

n - y:e-^T-)d = (1 - e^){(T;) - e^(r-,)} (6) 

and 



d 



(T,-) = (<^o|T^|0o) (7) 

Where r is the Bragg wave- vector, is the Dy 4f ground state wavefunction. 
Equations (5) to (7) show that A\{J, M) and A\{J^ M) are proportional to the dipole 
moment {Tl), octupole moment (Tf), and the triakontadipole moment (Tf) while 
A\{J,M) is proportional to the quadrupole moment (Tj), the hexadecapole moment 
(T2), and the hexacontatetrapole moment (Tl) . 

The total intensity is represented by scattering from the dipole {T^) and 
quadrupole (T|) of the 4f shell and reflects the same azimuthal angle dependence 
observed in the dipole transition at the L3 edge [23 • The dependence on the 
higher rank 4f multipoles arises from the splitting of the core hole states, which 
results in the M dependence of the amplitudes of the different harmonic oscillators. 
Correspondingly, these amplitudes are influenced by contributions from higher rank 4f 
multipoles which would otherwise cancel due to symmetry imposed by the 3j symbols 
in (5). Therefore, the core hole splitting allows the measurement of higher rank 
multipoles through its spectral shape than would normally be expected for an El 
resonance. 



High-order Dy multipoles in antiferromagnetic -D2/-B2C2 



6 



Figure 2. The best fit to the energy dependence for (a) cr, and (b) n incident 
radiation in the AFM+AFQ phase of DyB2C2 at 12K 

4. Discussion 

The theory detailed in section 3 is used to fit the measured energy dependence for a 
and TT incident polarization via a nonlinear least squares optimization process. The 
best fit is shown in Figure [5] with the oscillator positions indicated with arrows. The 
absorption correction is applied to the calculated spectra. In order to constrain certain 
features of the fitting procedure, Fj = a^Fj where Og is a scaling parameter. 
The parameters obtained in the fit are shown in Table [T] Aj typically depicts the 
onset of an absorption edge and the value of = 0.22 calculated for the free ion 
[24] gives an excellent fit. We find that the quadrupole interaction Q5/2 — —0.42 eV, 
and the internal magnetic field Hint = —0.98 eV. Hint is largely independent of the 
choice of wavef unctions and linewidths, and its order of magnitude is in line with 
that expected from atomic calculations of 3(i-4/ exchange integrals [23] . Variations in 
(T^) can be partially compensated by the line widths to yield a similar result. The 
linewidth for time-even contributions Fj increases progressively from 0.7 eV for 

F| to 4.5 eV for F| ,3, and is approximately a factor of 2.8 larger for time-odd 

2 '1 2 1 2 '1 2 

contributions F j vj/^ jg normalized to 1 and the best fit for the values of gives 
4-4 = 0.22, = 0.07, = -0.34, *f = 0.77, and *f = -0.78. We find that 
the octupole and triakontadipole moments are the dominant time-odd contributions, 
rather than the dipole moment. This is in accord with observations that the elastic 
properties of DyB2C2 are similar to the related compound H0B2C2 for which a strong 
influence from octupolar order is inferred [14j . We note in particular this is the first 
direct measurement of a triakontadipole moment. 

It should be noted that the analysis in terms of a core-valence interaction is 
consistent with the multiplet structure interpretation. The M4 and M5 absorption 
structure is « lOeP^ and w 5eV wide respectively, which is mainly due to the strong Sd- 
4/ and 4/-4/ Coulomb and exchange interactions [24] . The strong electron correlation 
leads to non-diagonal matrix elements in the quantum numbers L, S, and J of the 
core hole. The levels of the same M connected by these non-diagonal elements excite 
partly coherently and this makes it difficult to separate the core hole states. 

We note that the empty states of the 4/ shell selected by the resonant diffraction 
process and their corresponding relative transition intensities are not necessarily the 
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Table 1. Fitted Parameters in the AFQ and AFM + AFQ phases of DyB2C2 

Parameter T = 18K (AFQ) T = 12K (AFM + AFQ) 



^5/2 


1 9S7 7 pT/ 




A3/2 


1324.0 


1323.9 eV 




1^ 




*t 


-0.2 


0.22 




0.3 


0.07 


iTfl 




— L).o4 


iTf3 




U. M 







-0.78 


Qj 


-0.41 eV 


-0.42 eV 


Hjnt 


— 


-0.98 ey 


rj 


0.22 


0.22 


i 5 1 1 

2 ' 1 2 


u.o ev 


u. ( e V 




2.7 eV 


2.8 eV 




5.3 


3.2 eF 


2 '1 2 


1.9 eV 


2.2 eV 


-•^3 13 

2 ' 1 2 


5.4 eF 


4.5 eF 






2 


-1-5 13 




7.9 eV 


-■-5 15 




9.1 el/ 


3 1 1 




6.2 eV 


1 1 3 




12.8 ey 



2 'I 2 I 

This value is fixed to 1 



same as for an absorption process. Resonant diffraction is sensitive to the difference 
between electronic states and the conditions for nonvanishing intensity depend only 
on local symmetries (of the resonant atom) while absorption spectroscopy is sensitive 
to the average difference between electron states, and the conditions for nonvanishing 
terms must reflect the global symmetry (Neumann's principle) [351 HHj. For example, 
the magnetic moment of an antiferromagnet is not visible at the absorption edge, but 
resonant diffraction may yield large intensities proportional to the magnetic moment 
[37] . In the present theoretical framework the core hole is decoupled from the 4/"+^ 
state and the effect of the 4/ multiplet structure is indirectly taken into account by 
the effective widths of the M distributions which exceed the intrinsic lifetime width 
and describe the data surprisingly well. In this way, much physical intuition regarding 
the associated atomic quantities of interest (i.e. the atomic tensors (T^)) is retained. 
This compares to recent numerical calculations |5S] which included a crystal field 
splitting of the 4/ states and the full multiplet structure, where a similar magnitude 
of Qj was found (although of different sign) but no information on the (T^) could be 
extracted from such a calculation. This work further demonstrates that including the 
intra-atomic coulomb and magnetic interactions is a viable approach to account for 
the multiple spectral features and their broad distribution in energy. 
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5. Conclusion 



DyB2C2 exhibits an ordered pliase witli botli antiferroniagnetic (AFM) and antiferro- 
orbital (AFQ) order below T-^ = 15.3 K. Soft resonant x-ray Bragg diffraction 
experiments were performed at the Dy M4^5 edges and the energy dependence of 
the space group forbidden(OOi) reflection was measured in this low temperature 
AFM+AFQ phase. The Dy 4f multipole moments of rank 1 (dipole) to 6 
(hexacontatetrapole) were measured and their magnitudes are determined to be 
between 7 and 78 % of the the quadrupole (rank 2) moment. The dominant time- 
odd contributions are of rank 3 and 5, indicating a strong influence from octupolar 
and triakontadipole multipoles. The energy dependence of the (OO5) reflection is 
modelled successfully by including an intra-atomic core-hole interaction parameterized 
by quadrupole and magnetic interactions. A pseudo-multiplet structure was also 
introduced, allowing for the intereference between oscillators that have different 
properties under time-reversal due to their having different excited-state lifetimes. 
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Appendix A. Explicit Relations for X, 

The quantity Xq appearing in Equation (1) is a spherical tensor constructed from 
components of the incident and reflected polarization vectors. It is described by 
Equation (66) in W, and the values appropriate to the current study are: 



2X1/2 

3 



X^ = sin(20), X2 = -i= cos(20), X| = 1 
/ 

TT a : 

A, = e , Ai = e 

^ 2 ' ^ 2 

a'lr : 

xl = -e'\ xl = --e' 

^ 2 ' 1 2 

Due to the inability to measure the polarization of the outgoing beam without severely 
attenuating the measured signal, the spectra in figures 1 and 2 are combinations of the 
rotated and unrotated channels. Specifically, = 1^' 



Appendix B. Dy site symmetry and magnetic structure 



The structure factor appropriate to the Dy site symmetry and magnetic ordering 

vI/^ = (l-e^){(r,-)-e^(T^,)} (B.l) 
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does not account for the slight rotations of the magnetic moments proposed in, for 
example, [29; . In the above case, we arrive at the following 

n = MT^y (B.2) 

= -2i{{T^)' - (Tf>") (B.3) 

where (. . .)' and (. . .)" represent the real and imaginary parts respectively. If we 
include the proposed rotations of ±9° away from the site symmetry, we arrive at the 
following structure factor 

VP- = (1 - e^e'^){{T^) - eWe^(Tf,)} (B.4) 

which gives 

*f = -Qmi{T^)" + ?,M{T^)" (B.5) 

■^l = -1.98i(rf)' + 1.69(rf )" (B.6) 

Appendix C. Reduced matrix elements R^{r,x) 

R^(r,x) is calculated for the Dy'^^ ground state of ^i?i5/2 using 

x(2a+l)(25+l)^(2y+l)| 



K r X 
b y 





I 






I 






b 






where the relation between W'-"''''' and W*^"^''-'^ is given by, 

is a unit double tensor in the defined in accord with Judd [21 , where a is the 
rank of the spin component and b is the rank of the orbital component. Its calculation 
is facilitated with the following equation: 



L+L+b+S+S+a+s+l 
a,L,S 



{raLS\W^''-''^\ra'L'S') =7j[L,L',b,S,S',a] ^ (-1) 

a,L,S 

^{ L L' Ijl^ S' 1} i'^"<^LS{\l'''^^){r-^^\}ra' L' S') (C.3) 

where {V''aLS{\l^~^aLS) is a coefficient of fractional parentage (CFP) and values 
for them are tabulated for almost all electron configurations of interest in Nielson & 
Koster [30], along with values of {e\\V{K)\\e') = {^{2S + 1)}'^^^ {dWW"''^ \\e'). s 
and I are the spin and orbital angular momentum of the equivalent electrons coupled 
together in this process, and in the case of 4/ electrons have the values 3 and ^ 
respectively. The results for the calculation of reduced matrix elements (r, x) for 
the Dy'^+ ion (which has = 5 holes) in the ground state are tabulated 

in tables [C2] and [Cll 
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Table CI. Reduced matrix elements R^{r,x) for ^Hi5 and nh = 5 . 



= 1, J = I 



r 









1 




2 




3 




x=0 







4 
21 




2 

15 












1 




15 









8 

525 \ 


/l7 
' 15 







2 







4 
21 




Tr 

35 







2 / 
49 V 


' 34 
105 


3 














2 , 
175 V 


/646 
195 







4 





















1938 
455 



K = 1,J=^ 

71 1 2 3 4 5 
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2 
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